Electronic structure and bond competition in the polar magnet PDVO3 
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Density functional electronic structure studies of tetragonal PDVO3 are reported. The results 
show a an important role for both Pb 6p - 2p and V d - p bonding, with an interplay between 
these. This is discussed in relation to the possibility of obtaining magnetoelectric behavior. 
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I. INTRODUCTION 

Magnetoelectric and specifically multiferroic magneto- 
electrics have attracted renewed attention, 1 ' 2 ' 3 4 although 
the field has a long historyi^iSiL2iS Nonetheless, bulk 
materials that provide strong coupling between electric 
and magnetic fields at room temperature remain elu- 
sive. Some of the difficulties in realizing such materi- 
als were reviewed by Hilli^> Perovskite and perovskite de- 
rived phases, ABO3, provide the main inorganic ferro- 
electric materials used in applications. These are based 
on a cubic lattice with A ions on the cube corners and 
B ions at the cube centers, inside corner sharing octa- 
hedral O cages, Ferroelectric perovskites can be classi- 
fied as A-site or £?-site driven according to which ion 
is primarily responsible for the polar instability. One 
strategy for synthesizing multiferroics is to substitute 
one of the perovskite sites with magnetic ions in ma- 
terials where the ferroelectricity is driven by the other 
site, as for example in BiFe03ii£ However, in both A 
and -B-site driven materials the ferroelectricity gener- 
ally involves cooperation between off-centerings of both 
cations, and furthermore covalency between unoccupied 
metal states and O p states is important for the ferro- 
electric instability>ii*i2iiiii2ii£ii& From this point of view, 
V 4+ may be a particularly interesting ion for multifer- 
roics. 

In oxides, V 4+ is often magnetic and, significantly, 
can occur either at the centers of octahedral cages, as 
in CaVC>3 and SrVC>3, in which case there are interesting 
metal insulator transitions of Mott character?iZii& or it 
can strongly bind to one O forming forming a vanadyl 
ion with very short bond length, ~ 1.55 - 1.60 A. These 
vanadyl containing compouds are local moment magnets 
due to the formation of spin 1/2 moments on the V ions, 
and are band insulators, The band gaps are formed by 
a combination of strong crystal field splitting associated 
with the short V-0 distance, and on-site Hund's cou- 
pling, which exchange splits the 1/2 filled V non-bonding 
orbitalii2iS£ Thus, perovskite vanadates offer the possi- 
bility of combining a magnetic B-site ion with insulating 
behavior and stereochemical activity on both the A and 
B-sites. This was recently realized by the high pressure 
synthesis of PbV0 3 M^. 

PbVC>3 is insulating and occurs in the tetragonal, 
non-centrosymmetric structure, spacegroup PA/mm, like 
PbTiC>3, but with a very large tetragonality, c/a ~ 
1.23. The structural distortion may be described as off- 
centering of Pb and V accompanied by a distortion of the 



VC>6 octahedral cage. This is in contrast to CaVC>3 and 
SrVC>3, which have ccntrosymmctric GdFeC>3 type struc- 
tures deriving from octahedral tilting. PbVC>3 shows no 
structural transitions from OK to the decomposition tem- 
perature. However, a tetragonal to cubic transition does 
occur under pressure, 22 indicating the possibility of ferro- 
electric behavior under pressure and perhaps by chemical 
modification. It is noteworthy that at 300K the short- 
est V-0 distance is 1.67 A, which is longer than in 
most other vanadyl containing oxides. This implies that 
the vanadyl bond is relatively weak in PbVC>3. Further- 
more, because of the connection between the short V-0 
bond and the moment formation, one may anticipate the 
possibility of coupling between magnetism and the struc- 
tural distortion, which in turn, if ferroelectric, would be 
coupled to external electric field. Significantly, the Pb 
displacement with respect to the center of its O cage 
is 0.92 A, and the V displacement is 0.60 A. These are 
large and of comparable magnitude, indicating both ro- 
bust ferroelectricity if switching could be achieved, and 
cooperativity between the two cation sites, which is a 
characteristic of good ferroelectric perovskites. 

Shpanchenko and co- workers^ reported local spin den- 
sity approximation (LSDA) electronic structure calcula- 
tions and an analysis of the bonding and magnetism. 
These calculations were based on the linearized muffin-tin 
orbital method with the atomic spheres approximation 
(LMTO-ASA). This method is well suited to magnetic 
transition metal oxides, but is suffers from shape and ba- 
sis set errors in open structures with highly non-cubic 
site symmetries, as is the case for PbV03. They also re- 
ported LSDA+U calculations for ferromagnetic ordering 
only, with the augmented planewave plus local orbital 
(APW+LO) method. The comparison of these two leads 
to the conclusion that the formation of the spin polarized 
vanadyl bond depends on Hubbard correlation effects. 
Here we report LDA calculations done with the gen- 
eral potential linearized augmented planewave (LAPW) 
method, which incorporates the full crystal potential 
and has a much better basis set than the LMTO-ASA 
method, especially for the interstitial regions^ We use 
these to discuss the formation of the electronic structure 
and the bonding of the compound. We find results that 
differ significantly from those of Ref. 0, with a different 
conclusion regarding the vanadyl bond. Our results com- 
pare well with very recent calculations of Uratani and co- 
workersiS* We find an interplay between Pb-0 and V-0 
covalency that may be of importance for magnetoelectric 
materials. 
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II. COMPUTATIONAL METHOD 

As mentioned, the calculations were done using the 
general potential LAPW method. Except as noted, the 
experimental 300K crystal structure reported by Belik 
and co-workersSS was used. The calculated LSDA forces 
on the atoms are small with this structure, supporting 
the experimental determination as well as the conclusion 
that the vanadyl bond is properly formed at the LSDA 
level. The largest force is on the apical O of 0.032 Ry/ao 
towards the V, with a smaller force of 0.017 Ry/ao draw- 
ing the V towards the apical O for the lowest energy (see 
below) antiferromagnetic ordering. The LAPW sphere 
radii were 2.00 do, 1.70 clq, and 1.35 do, respectively for 
Pb, V and O. Local orbitals^ were used to treat the 
high lying valence states of Pb and V, as well as to re- 
lax linearization errors. The basis sets consisted of ap- 
proximately 1100 LAPW functions plus local orbitals for 
the primitive cells and correspondingly larger sizes for 
supercells. Zone sampling during the iterations to self- 
consistency was done using the special points method, 
with 40 k-points in the irreducible wedge for the prim- 
itive cell. Tests with other choices were made. These 
showed that the calculations were well converged. Mag- 
netic ordering was studied using \/2 x \/2 x 2 supercells. 
In all cases, the core states were treated self-consistently 
and fully relativisically, while the valence states were cal- 
culated in a scalar relativistic approximation. 



III. MOMENT FORMATION 

Stable spin moments of 1 ps were found for all mag- 
netic configurations, consistent with local moment mag- 
netism. This differs from the calculations reported in Rcf. 
I2H where the F ordering showed only a small moment 
in contrast to the C ordering. As mentioned, they also 
did LSDA+U calculations for the ferromagnetic case, but 
with the APW+LO method instead of the LMTO-ASA 
method that they used for the other calculations in their 
paper. The implication of their results is that the Hub- 
bard U is essential for a stable spin polarized vanadyl 
bond. We find a stable spin polarized bond in accu- 
rate LSDA calculations, with the implication that the key 
physics is hybridization with O and on-site Hund's cou- 
pling. Calculations were done for an artificial non-spin- 
polarized CBjSG j MS well as ferromagnetic (F), antiferro- 
magnetic in the ab plane (C), antiferromagnetic along all 
three directions (G), and ferromagnetic ab planes stacked 
antiferromagnetically along the c-direction (A). The cal- 
culated energies are given in Table [I] 

We find the lowest energy state to be antiferromagnetic 
type C, qualitatively in agreement with the calculations 
of Ref. |2ll and with energies in good agreement with 
those of Ref. 0- We note that the energies are not 
well described by a nearest neighbor Heisenberg model. 
Both the C and G orderings have antiferromagnetic ab 
planes but opposite stacking along c. Similarly, F and 



TABLE I: Energetics of various spin configurations of PDVO3 
relative to the non-spin-polarized case. Energies are given in 
meV per formula unit. 



F 


G 


A 


C 


-111.2 


-127.3 


-91.8 


-127.8 



A have ferromagnetic ab planes, again differing in the c 
stacking. However, the energies of C and G differ by less 
than 1 meV per V, while those of F and A differ by ap- 
proximately 19 meV. Thus interactions at least to second 
neighbor are needed to reconcile the energies. Consider- 
ing this, as well as the extremely small energy difference 
between the C and G orderings, the actual ground state 
could be complex. One may speculate that this is the rea- 
son for non-observation of antiferromagnetic superlattice 
peaks in neutron diffraction. We note that there is ev- 
idence for a transition at ~120K seen both in thermal 
expansion^ and resistivity^!. 



IV. BAND STRUCTURE 

We now turn to the electronic structure, starting with 
the F ordering. The band structure is shown in Fig. 
^ and the corresponding density of states (DOS) and 
projections in Fig. The band structure, in order of in- 
creasing energy, consists of a split-off Pb 6s derived peak, 
approximately 9 eV below the Fermi energy followed by 
a manifold of O 2p derived bands, and then, following an- 
other gap, a manifold of metal derived bands. Initially, 
these are V d derived, but at higher energy (~ 3 eV) 
these change to include Pb p and other characters. 

The crystal field may be regarded as a combination 
of an octahedral crystal field, which splits the d levels 
into a lower lying t2 g group and a higher lying e g group. 
These are then further split by the off-centering. The 
main contribution to the crystal field in oxides is from 
metal - oxygen hybridization. Since the O 2p bands are 
below the V d bands, the V d orbitals are formally anti- 
bonding, so the the stronger the hybridization of a given 
d orbital with O, the higher it will be pushed. The result 
is a splitting of the e g manifold into two and a splitting of 
the t2 g manifold into a low lying d xy non-bonding band 
and a higher two-fold degenerate d xz and d yz derived 
manifold. This explains the multi-peak structure of the 
d band DOS. In the case of V 4+ , there is one d electron 
in the d xy orbital, which forms local moments due to the 
Hund's coupling. The exchange splitting for F ordering 
is 1.43 eV, and not sensitive to the magnetic ordering 
(see Fig. |3J). This is sufficient to fully polarize the d 
orbitals. For F ordering this yields a half-metallic state, 
spin moment 1 /ib- 

With the experimental structure of PbVOs and F or- 
dering, the crystal field splitting is almost, but not quite, 
sufficient to open a complete gap between the majority 
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FIG. 1: Band structure of ferromagnetic ordered PbVC>3. 
The top panel shows the majority spin bands and the bot- 
tom shows the minority spin. 



spin d xy band and the higher lying d bands. This is be- 
cause of the bandwidth of the d xy band. As may be seen 
in majority spin panel of Fig. ^ the d xy band has almost 
no dispersion along the T — Z line and has weak disper- 
sion along M — A implying that the main hopping is the 
ab plane. Following the dispersion along T — M , the d xy 
band maximum is 0.3 eV above the bottom of the other 
d bands. 

We did a structural relaxation for the F ordering. The 
resulting LSDA positional parameters are zv=0.5689, 
z O i=0.2158, and z O 2=0.6872, yielding a V - 01 bond 
length of 1.67 A, which is slightly less than the exper- 
imental value. The corresponding experimental values 
are z v =0.5668, z O i=0.2102, and z O 2=0.6889. The calcu- 
lated full symmetry a g Raman frequencies are 190 cm -1 
(Pb motion against the other atoms in the cell), 408 cm -1 
(cations vs. O), and 838 cm~ 3 (V - 01 bond stretching, 
reflecting the short bond length and stiff interaction). 

The effect of in-plane antiferromagnetism is to narrow 
the d xy band by disrupting nearest neighbor hopping. 
This is because the exchange spitting is large. The result 
is a band gap as shown in Fig. [31 Because the crystal field 
splitting is sensitive to the exact V - apical O distance, 




FIG. 2: (color online) Electronic DOS and projections onto 
LAPW spheres for ferromagnetic PbVOa. The top panel 
shows the total DOS and projections onto the V d charac- 
ter, while the bottom shows Pb contributions. The middle 
panel is a blow up around Ef- Note that since the Pb or- 
bitals are more extended than the 2ao sphere radii, the Pb 
projection is proportional to the corresponding contributions, 
but reduced in magnitude. Majority and minority spin are 
shown above and below the axes, respectively. 
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FIG. 3: (color online) Electronic DOS of PbV0 3 for the lowest 
energy C-type antiferromagnetic ordering. 



the band gap should also be sensitive to this structural 
parameter. We can also understand in these terms why 
the energy difference between the G and C type orderings 
is small and that between the F and A type orderings is 
not, even though both pairs differ only in the c direction 
magnetic stacking. Specifically, in the F and A orderings, 
the majority d xy manifold overlaps the other d states so 
that the moment is not of pure d xy character, and there 
are non-trivial dispersions in the c direction in the region 
of the zone where there is overlap, such as M — A of 
Fig. n These dispersions provide interaction along c, 
and because they consist of is metallic d band overlap in 
only one spin channel, the resulting spin interaction is 
strongly ferromagnetic. Thus F is strongly favored over 
A. On the other hand, with antiferromagnetic in-plane 
ordering, the d bands are narrower, and so the d xy band 
does not overlap the other d states. Thus the c axis 
exchange coupling is much weaker. At the experimental 
structure it is very weakly ferromagnetic, favoring C over 
G, but if the gap were larger, either due to a shortening of 
the V - apical O distance or on-site Coulomb repulsions, 
one would expect it to become antiferromagnetic. In this 
case, if the magnetic structure is simple, the ground state 
would become G type. Thus, at least at the LSDA level, 
PbV03 is near a structure dependent cross-over between 
two magnetic states. 



V. ROLE OF PB 

Turning now to the role of Pb, the bottom panel of 
Fig. |2]shows projections of the DOS onto the Pb LAPW 
sphere (radius 2 ao). The large peak below the O 2p 
bands is from the Pb 6s states. As may be seen, these mix 
with the occupied O 2p bands, but play almost no role 
in the unoccupied conduction bands. Mixing of occupied 
states does not produce bonding. On the other hand, the 
Pb 6p states are low lying in the conduction bands and 




FIG. 4: (color online) Electronic DOS of hypothetical ferro- 
magnetic CaVOa constrained to the PbV03 crystal structure 
(see text) and a blow-up around Ef- 



mix with the O 2p bands, which are the same bands that 
hybridize with V d states. Thus there is a mechanism for 
competition between the vanadyl bonding and the Pb - 
O bonding. 

We emphasize that while the lone 6s pair of Pb 2+ is 
often discussed in terms of the stereochemical activity 
of Pb, these orbitals are not responsible for the activity 
directly. Rather it is that in T1(I), Pb(II) and Bi(III) 
the atomic structure that leads to the position of the 6s 
orbital also leads to 6p levels that are low lying in the 
conduction bands and very spatially extended favoring 
hybridization with occupied ligand orbitals. We empha- 
size that while the lone 6s pair of Pb 2+ is often discussed 
in terms of the stereochemical activity of Pb, these or- 
bitals are not responsible for this activity directly. Rather 
it is that in Pb(II) and Bi(III) the atomic structure that 
leads to the position of the 6s orbital also leads to 6p lev- 
els that are low lying in the conduction bands and very 
spatially extended favoring hybridization with occupied 
ligand orbitals. 

Sr 2+ is similar in size to Pb 2+ , but does not have the 
lone pair and the resulting low lying spatially extended 
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unoccupied p states. As a result, it would not form struc- 
tures with short bond lengths to O like the reported 
structure of PbVC>3. Ca 2+ is a smaller ion, also lack- 
ing the lone pair, and is more compatible with the bond 
lengths in PbVC>3, though, as mentioned, CaVC>3 actu- 
ally forms a centrosymmetric tilted GdFeC>3 type struc- 
ture. As such, we did calculations for hypothetical ferro- 
magnetic ordered CaV03 with the atomic positions fixed 
at those of PbVC>3 in order to elucidate the role of Pb. 
The calculated DOS is shown in Fig. As may be seen, 
the pseudogap between the majority spin d xy band and 
the remaining d bands is reduced and the half-metallic 
character is lost (the spin moment is 0.88 hb instead of 
1 He)- As mentioned, in PbVC>3 the calculated forces on 
the atoms are small. Keeping the structure fixed, there is 
a strong force of 0.08 Ry/cio pushing the V away from the 
apical O, i.e. disrupting the vanadyl bond. Thus indeed 
there is substantial coupling between the Pb-0 interac- 
tion and the V-0 interaction, seen both in the stability 
of the vanadyl bond and the magnetic properties. This 
is due to bond competition with O. Specifically, the re- 
moval of Pb leads to broadening of the ti g bands, which 
which spoils the crystal field scheme that stabilizes the 



vanadyl bond. 



VI. SUMMARY 

To summarize, LSDA calculations support the exper- 
imental conjecture that PbVC>3 may be the basis for in- 
teresting magnetoelectric materials. This is based on the 
electronic coupling between Pb and V due to mutual 
bonding with O and the stabilization of magnetism by 
a mechanism connected with the V-0 bonding. The 
key challenge is to find chemical modifications that make 
the material more switchable, so that ferroelectricity can 
be observed. 
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